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ABSTRACT 


The  theorem  of  least  work  was  used  to  investigate 
the  stress  pattern  and  bowing,  i.e.  uniform  curvature, 
of  an  initially  straight  thin  tube  subjected  to  a  steady 
temperature  distribution  independent  of  the  longitudinal 
coordinate  of  the  tube.  Simple  equations  were  developed 
for  stress,  strain,  and  bowing.  They  show  that  if  the 
thermal  strain  is  expressed  as  a  Fourier  series^  the  tan¬ 
gential  stress  is  small  and  is  a  function  of  the  first 
harmonic  only,  the  axial  stress  is  a  function  only  of 
the  second  and  higher  harmonics,  and  the  bowing  is  a 
function  only  of  the  first  harmonic . 
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1.  Nomenclature. 

A  =  cross  section  area  of  tube  wall 

a,b  =  coefficients  used  to  specify  axial  strain 
(See  Equation  4) 

E(e)  =  local  thermal  strain,  presumed  a  function  of  Q 
only,  not  of  r  or  z. 

coefficients  in  Fourier  expansion  of  £(&) 

E  =  Young's  modulus  of  elasticity 

n  =  an  integer  denoting  the  nth  term  of  a 
Fourier  series 

r  =  radial  coordinate,  measured  positive  outward 
from  mid-thickness  of  tube 

ij,  =  tube  radius  to  mid-thickness  of  wall 

t  =  tube  wall  thickness 

T(6)  =  local  temperature 

=  base  temperature  from  which  thermal  strain  is 
measured 

U  =  strain  energy  per  unit  length  of  tube 

u(@)  =  local  radial  displacement  of  mid-thickness 
curve  of  tube 

u„,un  =  coefficients  in  Fourier  expansion  of  u(@) 

v(@)  =  local  tangential  displacement  of  mid-thickness 
curve  of  tube 

vn,vM  =  coefficients  of  Fourier  expansion  of  v(@) 
x,y,z  =  tube  coordinates  (see  Figure  1.) 

**  =  coefficient  of  linear  thermal  expansion 

ft  =  ^  + 

~r6  =  elastic  strain  components  in  tangential  and 

longitudinal  directions,  respectively 
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j  =  total  strain  (elastic  plus  thermal)  in  the 
tangential  and  longitudinal  directions 


respectively 

t  + 


n 

9  = 


angular  coordinate  measured  clockwise  from 
the  positive  z-axis,  (see  Figure  1). 


angular  position  of  the  plane  of  bowing 
K  =  curvature  of  bowed  tube 
v  =  Poisson's  ratio 


£=/--£- 
<  !Z  c 


P,  =  / 


/2, 


direct  stress  in  the  radial,  tangential,  and 
longitudinal  directions,  respectively 


shear  stress  in  the  r#~plane 

i+vCi-ftt ) 


Tr& 

#,  =  #  P 

X 

t 

<p, 


rj 

T  ^  \  E 


 /- 


/- 


Cu  = 


_  d  tv 

cfe 


TJ 


where  ou  is  any  representative  variable 


2.  Introduction 

The  problem  of  stress  distribution  in  a  tube  with 
temperature  variation  independent  of  the  axial  coordinate 
has  been  solved  under  the  following  two  sets  of  circum¬ 
stances:  (1)  an  exact  solution  [2]*  in  the  sense  of  the 
theory  of  elasticity  for  steady-state  temperature  config¬ 
uration,  and  (2)  a  very  elementary  soiution[5j  in  the 
framework  of  elementary  strength  of  materials .  Solution 
(1)  is  of  limited  usefulness,  not  only  because  of  its 
limitation  to  the  steady-state  configuration  but  also 
because  of  its  mathematical  complexity.,  Solution  (2)  was 
developed  to  investigate  technically  important  cases  of 
the  bowing  of  pipes  partially  filled  with  a  cryogenic 
fluid,  such  as  liquid  oxygen.  This  solution  assumes 
that  only  longitudinal  strains  and  stresses  in  the  pipe 
are  an  important  cause  of  bowing .  This  in  effect  is 
taking  Poisson’s  ratio  equal  to  zero.  Presuming  that 
Poisson’s  ratio  is  zero,  while  actually  it  is  in  the 
neighborhood  of  0.3  for  most  piping  materials,  led  to  a 
simple  and  useful  solution  for  bowing  radius  and  axial 
stress  but  said  nothing  about  tangential  stress.  The 
degree  of  approximation  involved  was  quite  unknown;  there¬ 
fore,  it  was  the  purpose  of  the  present  investigation  to 
obtain  a  more  accurate  solution  than  that  of  Flieder, 
Loria,  and  Snith  [53  without  involving  the  mathematical  com- 

^Numbers  in  brackets  refer  to  works  listed  in  the 
bibliography,  page  45* 


3 


Then  strain  energy  is  differentiated  partially  with 
respect  to  each  such  coefficient,  each  such  partial  deri¬ 
vative  is  equated  to  zero,  and  finally  the  resulting  set 
of  equations  is  solved  for  the  corresponding  values  of 
the  coefficients. 

This  procedure,  which  was  used  for  the  present  investi¬ 
gation,  gives  the  displacement  which  minimizes  the  strain 
energy  (subject  to  the  implied  constraints  of  the  assumed 
displacement)  and  thus  corresponds  to  equilibrium. 

That  the  theorem  of  least  work  can  be  applied  to  the 
problem  at  hand,  where  temperature  changes  are  involved, 
is  not  obvious.  Argyris[l]  and  Boley  and  Weiner  [2]  have 
discussed,  in  general  terms,  the  energy  theorems  and  thermal 
stresses  where  temperature  changes  and  non-linear  elasti¬ 
city  apply.  Even  so  it  seems  appropriate  to  establish, 
by  using  the  following  physical  considerations,  that  the 
theory  of  least  work  can  be  applied  to  the  present  problem. 

Consider  the  structure  cut  into  a  number  of  elemental 
volumes  such  that  the  temperature  change  is  uniform  in 
each  element.  A  change  in  temperature  causes  thermal 
strains  but  not  stresses;  however,  the  elements  no  longer 
fit  together  continuously.  In  order  to  fit  the  elements 
together  into  a  continuous  structure  having  an  assumed  dis¬ 
placement  pattern  from  the  original  shape,  stresses  must 
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be  applied  which  cause  elastic  strains  in  accordance  with 
the  theory  of  elasticity.  Since  the  stresses  are  applied 
holding  temperature  constant,  the  strain  energy  can  be 
calculated  using  the  equation  previously  cited,  f<rd~r 
The  values  of  the  elastic  strains  are  determined  by  the 
difference  between  the  total  strain,  determined  from 
geometric  considerations  of  the  assumed  displacement  pat¬ 
tern,  and  the  thermal  strains  resulting  from  known  tempera¬ 
ture  changes.  The  theory  of  least  work  now  applies  and 
states  that  the  strain  energy  is  a  minimum  when  the  assumed 
displacement  pattern  is  that  corresponding  to  equilibrium. 
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3.  Analysis  of  Stress  and  Strain. 

In  order  to  proceed  with  the  analysis,  certain  condi¬ 
tions  must  be  met  and  assumptions  made,,  The  tube  must  be 
constructed  of  an  isotropic  elastic  material  and  the 
temperature  distribution  must  be  a  known  function  of  the 
angular  coordinate,  6  ,  and  independent  of  time  and  the 
longitudinal  coordinate,  z.  This  is  the  resulting  tempera¬ 
ture  distribution  when  a  tube  is  partially  filled  with  a 
cold  liquid.  Also  no  external  forces  or  moments  may  be 
applied  to  the  tube.  The  necessary  assumptions  are  as 
follows: 

1.  The  tube  is  ,,thinn  and  thus  and  Trg  are  small 
and  can  safely  be  neglected  in  the  analysis. 

2.  The  temperature  is  constant  across  the  tube  wall 
thickness. 

3.  Young’s  modulus  of  elasticity  and  Poisson’s  ratio 
are  constant.  This  is  a  reasonable  assumption  as  shown  by 
values  for  302  stainless  steel  in  Figures  1  and  2. 

Stainless  steel  is  typical  of  most  metals. 

4.  Plane  normal  cross  sections  remain  plane. 


By  neglecting  <Tr and  equilibrium  cannot  be  satisfied 
and  any  results  are  necessarily  approximate.  A  similar 
assumption  has  been  made  by  others  under  similar  circumstances 
and  the  solutions  obtained  were  acceptable  when  compared 
with  the  exact  solutions  according  to  the  theory  of  elasti¬ 
city.  See  DenHartog[4J  pages  221  to  223  and  page  240.  The 
latter  portion  of  Appendix  III  shows  the  relative  magnitudes 
which  may  be  expected  of  <rr  and  rrg  . 
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5»  In  accordance  with  Saint  Venant’s  principle  the 
stress  pattern  exists  uniformly  along  the  tube  except  at 
short  distances  from  the  ends. 


-400  -300  - ZOO  -too  O  /OO  200  300  400 


TEMPER ATURE/f 


Figure  1 „  Modulus  of  Elasticity  of  AISI  302 
Stainless  Steel 


Figure  2.  Poisson4 s  Ratio  for  AISI  302 
Stainless  Steel 


From  T.  F.  Durham,  R„  M.  McClintock,  and  R„  P„ 
Reed,  ’’Cryogenic  Materials  Data  Handbook", 
National  Bureau  of  Standards,  PB  Report  171809, 
1961. 
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A.  Symmetric  Case. 

Temporarily  the  temperature  distribution  will  be 
limited  to  one  symmetric  about  the  y-axis.  The  coordinate 
system  used  is  shown  in  Figure  3. 


x 


Figure  3 »  Typical  Section  of  Tube 


When  there  is  a  change  in  temperature  of  the  tube,  a 
point  P  located  on  the  mid-thickness  curve  will,  in  gen¬ 
eral,  be  displaced  to  a  new  location  Pi  The  radial  dis¬ 
placement,  u,  and  the  tangential  displacement,  v,  are 
positive  in  the  directions  indicated  in  Figure  4° 


Figure  4»  Displacement  of  Mid-thickness  Curve 
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The  tangential  strain  at  any  point  can  be  expressed  as  the 
strain  of  the  mid-thickness  curve  plus  the  strain  due  to 
bending  of  the  tube  wall,  which  in  turn  can  be  expressed 
in  terms  of  the  radial  and  tangential  displacements  as 

9  r  ~  r&~  -  (1) 


follows 

fo 


For  derivation  see  Appendix  I, 

Since  the  temperature  distribution  is  symmetric  about 
0=0,  the  radial  displacement  will  be  an  even  function  of 
$  and  the  tangential  displacement  will  be  an  odd  function 
of  0  .  Expressed  as  Fourier  series  they  are 

op 

u-  Y  cos(n&) 


As  o 
oo 


(2a) 


(2b) 


and  (n9'>  . 

H  =  / 

The  necessary  derivatives  are 

u  =  -Z  nZun  cos(ne) 

V -  Y  n  ^  cos(n&)  . 

*  -  j 

Substituting  Equations  2  into  Equation  1  we  obtain 

j:  !  +  +  n  IK,]  cos  (*  &)  .  (3) 

The  assumption  that  plane  cross  sections  remain  plane 
requires  that  the  axial  strain  be 


=  a  +  h(r0  +  r)  cosfe)  ) 


(4) 


where  a  and  b  are  constants  to  be  determined, 
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Because  there  are  no  external  forces  or  moments  applied 


to  the  tube  dA  —  0  (5) 

£y<r*dA=0  (6) 

fA*<*idA=0  (7) 

where  y  =  C r*  +  r) cos  6  and  d A  *  +  r)  dr  d  Q  .  (g) 


Equation  7  is  automatically  satisfied  as  a  consequence 
of  symmetry. 

With  the  assumption  that  radial  stress  is  zero, the 
stress  condition  in  the  tube  is  that  of  plane  stress. 

The  equations  of  plane  stress  are 

<re=  [et+ve'-O+^SM]  (9) 

yf?de*  +  v€*~°*‘')£<'e)]  ,  <io) 


where  £(&)is  the  thermal  strain.  Expressed  as  a  Fourier 

OQ> 

series  g(e)  =  8a  +  2  c°s  (n9)  ,  (11) 


where  S, 


Z TT 


iir  ,.  T 

©e  (T)  d  T 

-O  "T, 


and 


j,ifr 

£*  *  j  <xCr)dT  cos (nB)dB  •  n~ 

o  T-g 


The  thermal  strain  is  an  even  function  as  a  consequence  of 

the  temperature  being  symmetric  about  the  y-axis.  In 
rT(«) 

general  E(d)~  J  c*(T)J~f  0  Using  the  above  expression 
tb 
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for  axial  stress  and  integrating  Equations  5  and  6  permits 
determining  a  and  b  in  terms  of  displacement : 

«  =  <7^K>£0-  vJji*  (12) 

Q 

b-fi[a +»)£,-  (13) 

where  £  =  /  ~  j^*.  ,/S*  /  +  and  £  =  /  +  ^  . 

The  integration  is  performed  in  Appendix  II . 

The  elastic  strain  energy  for  the  case  of  plane  stress 
as  discussed  on  page  4  is 

u  =  if{L€*-£te)]ok  +  [ee-£(6)]<fe}  c/ >4  (14) 

per  unit  length  of  tube. 

Substituting  Equations  9  and  10  for  stress  into  Equation 
4,  multiplying  and  collecting  terms,  we  get 

=  Jig+et-nve^-z-O+vlECSiese,]  +zO+»)£(#f]jA.  U5> 

Integrating  Equation  15,  we  get  the  following  expression: 

[z?  ^ + n2  J 

+  i*r‘[2f4(io+/?b(“i+'/t)]  iE0+fi&  b£t] 

-  4f?  &)£*]  +zO+v)\z£.+'Ls5  ]'  (16) 

See  Appendix  II  for  details  of  the  integration. 
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Applying  the  theorem  of  least  work  as  discussed  pre¬ 
viously,  we  are  now  in  a  position  to  evaluate  the  coeffi¬ 
cients  uft  and  vn  by  taking  the  partial  derivative  of  'X 
with  respect  to  each  coefficient  and  equating  the  result 
to  zero. 

For  the  coefficient  u„, 


~4(i+i>)£jc-%£ )  ^  O  „ 


Therefore,  ^  (//or0  S0 


where  =  — — £ -  . 

r*  /->*<$ 


(17) 

(Id) 


For  coefficients  u,  and  v,  , 

X-  r0£bZ  +  -£*.&,+/,)*  +  &J?pb(u,+  v;)  -Z6+1/)/6k  b£t 

-zG+v)  (u,+  i/j£, 
c 


where  only  those  terir^s  that  contain  uf  or  v7  have  been 
written.,  Note  that  b  contains  the  sum  (  u7  +v#  ) .  It  is 

not  possible  to  evaluate  uf  and  v,  separately  because 


equals  c)X  0  However,  it  is  sufficient  to  obtain  their 

^  a  V, 

sum  because  that  is  the  only  form  in  which  u,  and  v,  appear 
as  can  be  seen  from  Equation  19 «  Taking  and  equating 
to  zero: 


+ zm.  +  -xnMt  [i-js?]*o. ,  20) 
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Solving  for  (u,+v,  )  we  get 

i 

(“,+  K)  -  P,  r0  S,  where  ^ 

For  all  other  coefficients,  u„  and  vn  ,  n>l  : 

jt  -  "  v«fa-0+»)f>«roah=  0% 


(21) 


and  +  nVn-0+v)ro£n= 


0 


yields  n  VA  ^  (l+P)fc6n  ~ /°n  &n  . 


(22) 

(23) 

(24) 


Substituting  Equation  24  into  22  yields  Cfn-fi)u^o  . 

But  C?n~f*  )  ^  0  ,  therefore, 

CL„~0  ;  VA=  0±vhs_£iL  t  (25a,  b) 

The  radial  and  tangential  displacements  are  now 
known  as  functions  of  the  thermal  strain,  and  thus  tempera¬ 
ture.  Substituting  Equations  IS  and  21  into  Equations  12 
and  13  respectively,  yields: 


cl- 


•—  where  £  —  1  +  vO  f  (pQ)  (26) 

O 

and  b  =  j*,  £j  where  “ -^-[f  +  v0~$)] m  (27) 


Thus  the  axial  strain  is 


e2  “  $0£.  +0+i )$,£,c°s(e). 


(28) 
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Substituting  Equations  18,  21,  and  25  into  3  gives  the 
tangential  strain, 

ao 

o  +  8, c*s  +  0+ k)X  6ncos(ne)  ,  ( 29 ) 

H-Z 

Remembering  that  £&> «  £0  + E.cosCe )  +J*  Encos(ne) 

*~2j 

0O_ 

then  Z  8„  cos  (nO)  =  8  (6)  ~50  S,  cos  8  and, 

ft-Z. 

et’‘Qtr')-v]£.+M-i)-v]e,cas6  +0*y)£(e).  (30) 

Substituting  these  expressions! or  strain  into  the  stress 
formulas,  Equations  9  and  10,  determines  the  stress: 

<5^=  -E£( a)+  +0-»l)l&>+l{$r,)+ v@r')+0-»,j+$likcos&}  ( 3 1 ) 

and  <%=^M&)+»($T$ao+l(</'t-/)+v($rO  +  »$tj:]E,eosGj%  (32) 
This  concludes  the  stress  analysis  for  the  symmetric  case,, 

B.  Anti -Symmetric  Case* 

If  the  temperature  distribution  is  anti “symmetric 
with  respect  to  0=0,  the  thermal  strain  is  anti-symmetric „ 
Expressed  as  a  Fourier  series 

OO 

6  Ce)  =Z  8n  stn  Cne)  (33) 
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where 


Note  that  there  is  no  constant  term  in  the  Fourier  series. 

It  is  seen  that  this  is  similar  to  the  symmetric  case 
rotated  90  degrees  and  we  can  apply  some  of  the  results  of 
the  previous  case. 

The  first  harmonic  is  the  only  one  contributing  to  axial 
strain,  so  the  axis  of  rotation  is  the  y-axis  and  the 
longitudinal  strain  can  be  written  as 


=  a,  +  b,  (r0-+rjs/n  6  >  (34) 

where  a, is  automatically  zero  because  of  anti-symmetry. 

The  radial  displacement  will  be  an  odd  function  and  the 
tangential  displacement  an  even  function  of  the  angular 
coordinate.  Expressed  as  Fourier  series  they  are 

OO 

a.(9)  =*£_  &„s/n  (no)  (35) 

/ 

oo 

and  V(9)=JE  V„  cos(nO)  .  (36) 

Note  that  v  contains  no  constant  term.  A  constant  term 
would  imply  a  rotation  of  the  tube. 

Substituting  Equations  35  and  26  into  3»  we  obtain 


€^=  -£j[[[i+£(n-0]u„  -  n  Vu'jsmCne)  .  (37) 
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Letting  ~  Vn  —  ^  n  yields 

\  +nXn}sm(nd)  .  OS) 

The  expressions  f  or  ^  and  are  of  the  same  form  as 
in  the  symmetric  case  except  that  sin(r?0)  is  substituted 
for  cos{n9) «  The  only  integrals  that  do  not  vanish  when 
integrated  over  the  interval  0  to  2tt  are  Js//j  Z(nQ)d &  <, 

These  integrals  have  the  same  values  as  fco,sx(nO)d@  for 
the  interval.  Thus  by  making  the  substitution  and 

writing  si n(n&)  for  cos(«0),  the  solution  of  the  symmetric 
case  is  applicable  and  the  results  are  as  follows: 


Q+n  )$,£,  sin  6  (39) 

ee=  Sin 9  +  (l+V)E (0)  (40) 

c rz=  -£E(d)  +  +v(ft -/)  +  £ ]st sin  d  (4D 

<r6  =  jfp:  fa,  ~0 4  v($rO  +»§,  -fo  ]if  s/n  o .  ( 42 ) 


The  stress  analysis  is  now  complete  for  the  anti-symmetric 
case . 
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C-  General  Case- 

For  a  general  temperature  distribution  the  thermal 
strain  can  be  expressed  as  the  sum  of  the  symmetric  and 


anti-symmetric  cases,  an  even  and  an  odd  function- 

°°  O© 

Thus  £(&)  =  €0 +2l£flcos(n6)  +  ]E  €nsm  (n&)  )  (43) 

a  complete  Fourier  series-  The  general  solution  is  as 
follows:  €.i=<&0£o+0+j:)$j£jcosO-i-Sls/n@]  (44) 

€e  =  +[(<£, -/) ->][£, cos 8  o]  +(/+») £(e)  (45) 

(%-~E  E(6)  +  ~jjr{t($o~0  +  »(farl)  +  0-^)]6o 

+  [£$,-0  4  vfy,  ~0*0-^  *$,f][£/cas  &  +£,sw0]}’  (46) 


The  stress  and  strain  resulting  from  a  general  temperature 
distribution  is  now  known- 

D-  Simplified  Equations  for  Stress  and  Strain- 

It  is  instructive  to  consider  the  magnitude  of  §t  , 

(po  ,  §t  ,  and  lpt  in  order  to  simplify  the  equations  of 
stress  and  strain-  The  stress  and  strain  equations  (Equa- 
tions  45  >  46  and  47)  contain  the  terms  {§0-/  ),  (^-/  ), 

(<^-/  ),  and  (^-/  )  which  vanish  if  the  ?s  and  ^?s  are 
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unity o  Table  I  shows  that  for  tubes  of  practical  wall 
thickness  this  is  very  nearly  so. 

TABLE  I 


The  Effect  of  Wall  Thickness  on  Stress  and  Strain 


0.1 

0.05 

0.01 

j". 

0.99992 

0.99996 

1.00000 

1. 00027 

1.00012 

1.00005 

io 

1.00027 

1.00014 

1.00001 

$, 

0.99825 

0.99913 

0.99972 

for  v  =  0.296 


Substituting  unity  for  ^  ,  $>f ,  and 


Equations  44 ,  45,  46,  and  47  become: 

€z  =  6o+0+-£ )(S,cos9  +€,S/ri  & )  (43) 

ed=  6+»)£(e)  -p(E*  +£,cos  e  +£,  s//'?  <9  J  (49) 

OJ  =  -e\e(Q) +  (50) 

0"  =^T ^  fey  cos  £  +  £, 9)  .  (51) 


These  equations  should  be  sufficiently  accurate  for  most 
design  applications . 
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4.  Bowing  of  Tube. 


A  tube  subjected  to  the  assumed  temperature  distribu¬ 
tion  will  bow  into  the  arc  of  a  circle,  or  if  restrained 
develop  stresses  equal  to  those  required  to  straighten  the 
bowed  tube.  This  is  of  practical  concern  in  pipe  design 
and  has  been  discussed  by  Flieder,  Smith,  and  Wetmore[6j. 
From  the  analysis  of  the  bending  of  an  initially  straight 
beam  [7]  it  is  known  that  the  curvature,  K  ,  is^-  ;  where 
c  is  the  distance  from  the  neutral  axis  and  €  is  the  strain 
resulting  from  bending.  The  axial  strain  on  the  mid- thick¬ 
ness  curve  of  the  tube  is  &+£,&/»  6)  ?  obtained 

by  taking  r=0  in  Equation  44 .  The  term  denotes  a 
general  elongation  of  the  tube  and  does  not  contribute  to 
bowing.  The  S,  c°£  &  term  results  in  bowing  in  the  yz-= 
plane.  The  strain,  €  ,  at  ®  *0  is  $>£>,  at  a  distance  r0 

from  the  neutral  axis.  Thus ,  Kut~  is  the  curvature 

in  the  yz-plane.  The  £>t6.,s/n9  term  results  in  bowing  in 
the  xz-plane.  The  strain  at  is  and  as  before  the 

distance  to  the  neutral  axis  is  ri  .  Therefore,  Kxts 


is  the  curvature  in  the  xz-plane.  The  curvature  in  general, 


(52) 


is  the  vector  sum  of  and  Ayz-,  and  the  plane  of  bowing 


is  given  by  dH  -  , 


(53) 
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measured  clockwise  from  the  y-axis.  Again  it  is  possible 
to  set  =1  without  significant  error  so  that 


r0 


(54) 
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5.  Results,  Design  Equations. 

The  analysis  leads  to  the  following  simplified  equa¬ 
tions  which  can  be  used  to  design  thin  tubing  systems 
subjected  to  steady  temperature  distributions  that  are  in¬ 
dependent  of  tube  length. 

1.  Axial  strain 


€a  =  +  )(E,cos&  +E,s/n&) 


(4*) 


2 .  Tangential  strain 


-v(eo+s,cos0 +£,,sme) 


(49) 


3.  Axial  stress 


4.  Tangential  stress 


ft  cos0  +  8,s/na) 


(51) 


5 .  Bowing  curvature 


(54)&(53) 


rvc  TC&J 


Where 


z7r  S  T($J 


£m  =  oc(r)dT  cos(nO)(Jg 


°  Ts 


_2X  T(Q) 


E«~Z7rf  f  <*(r)dT  sm(n&)  dd 


°  tb 
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6.  Discussion  of  Results. 

When  thermal  strain  is  expanded  in  a  Fourier  series, 
the  tangential  stress  and  bowing  as  given  by  Equations  51 
and  54  are  functions  of  the  first  harmonic  only,  and  the 
axial  stress  as  given  by  Equation  50  is  a  function  of  the 
second  and  higher  harmonics.  To  interpret  the  physical 
meaning  of  this,  consider  the  case  where  £(o) = £, 

When  the  thermal  strain  is  projected  onto  a  yz-plane,the 
first  harmonic  transforms  to  a  linear  function  and  the 
second  harmonic  transforms  to  an  even  function.  It  is 
known  [8]  that  linear  thermal  strain  does  not  cause  stress 
in  a  free  beam,  but  causes  bowing,  and  that  an  even  thermal 
strain  function  causes  stress  but  not  bowing.  Thus,  the 
results  are  those  expected  from  superimposing  elementary 
solutions  for  simple  beams. 

Equations  44  thru  47  include  a  first  approximation  for 
the  effect  of  wall  thickness;  however,  the  basic  assumptions 
of  this  analysis,  particularly  that  the  temperature  does 
not  vary  through  the  wall,  suggest  that  the  results  be 
confined  to  thin  walled  tubes.  Accordingly,  Equations  4& 
thru  51  and  54  are  suggested  as  being  the  appropriate  ones 
to  use  for  actual  design  and  analysis.  These  results 
afford  a  slight  extension  of  those  obtained  by  Flieder, 
Loria,  and  Smith  in  the  following  respects:  (1)  they  give 
a  variation  of  axial  stress  through  the  wall  thickness  and 
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(2)  they  give  a  value  for  tangential  stress. 

The  present  results  and  those  of  Flieder,  Loria,  and 
Smith  are  comparable  for  the  case  where  E,  V  ,  and  oc  are 
constant.  When  these  conditions  are  imposed,  Equation  50 
is  identical  at  the  wall  mid-thickness  to  the  result  devel¬ 
oped  by  Flieder,  Loria,  and  Smith.  Note  that  0^  is 
independent  of  Poisson’s  ratio.  Equation  50  gives  a  minor 
variation  of  through  the  wall;  they  did  not  consider 
any  radial  variation  of  .  Flieder,  Loria,  and  Smith  do 
not  give  a  result  for  tangential  stress.  Equation  51  shows 
that  tangential  stress  is  zero  at  the  tube  wall  mid-thickness 
and  that  it  is  small  for  tubes  with  slight  bowing.  Equation 
51  also  shows  that  tangential  stress  vanishes  if  Poisson’s 
ratio  is  zero. 

Comparison  of  the  results  for  bowing  show  that  there 
is  a  very  slight  effect  due  to  finite  wall  thickness  and 
non-zero  Poisson’s  ratio  as  given  by  Equation  52.  However, 
the  simplified  equation,  Equation  54»  which  is  here  recom¬ 
mended  for  purposes  of  design  and  analysis,  is  essentially 
identical  to  that  given  by  Flieder,  Loria,  and  Smith. 

Essentially,  we  have  shown  that  influences  neglected 
in  the  analysis  of  Flieder,  Loria,  and  Smith  have  only  the 
very  slightest  effect  on  the  final  result.  However,  the 
expression  of  the  thermal  strain  as  a  Fourier  series  has 
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simplified  the  final  form  of  the  equations  and  aids  in  the 
visualization  of  the  effect  that  a  given  temperature  dis¬ 
tribution  may  be  expected  to  have  on  a  tube. 
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simplified  the  final  form  of  the  equations  and  aids  in  the 
visualization  of  the  effect  that  a  given  temperature  dis- 
tribution  may  be  expected  to  have  on  a  tube. 
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APPENDIX  I 


Derivation  of  Expression  for  Tangential  Strain 

The  temperature  change  from  point  to  point  around  the 
circumference  of  the  tube  results  in  thermal  strains.  There 
is  assumed  to  be  no  temperature  gradient  across  the  wall 
thickness  of  the  tube,  thus  the  thermal  strain  is  indepen¬ 
dent  of  radius.  Superimposed  upon  this  strain  is  another 
strain  which  varies  with  radius  and  is  caused  by  the  change 
in  curvature  of  the  tube  wall. 

To  derive  an  expression  for  that  portion  of  the  strain 
independent  of  radius,  consider  two  points  A  and  P,  close 
together,  on  the  mid-thickness  curve.  When  the  temperature 
changes,  the  points  will  move  to  new  locations  A,  and  P(  as 
shown  in  Figure  5. 


x 


Figure  5.  Tangential  Strain  at  Mid-thickness. 
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Let  /•  and  7  be  unit  vectors  in  the  i  and  j  directions 


respectively.  The  bar  denotes  a  vector. 

Then  OA,  =  vl  +  (rQ+u)j  j  (la-1) 

oq  =  (r0+tz+do)\sm(dQ)l  +  cos(d&)j]  >  (la-2) 

QP,  =  (v+d v)[cos (c/0) 7  +  s/n(d6)j]  ?  ( la_3 ) 

and  A,f)  ~  OQ  QFf  ~~  OAf  t  (la-4) 

Because  dd  is  a  small  angle,  sin (dd)  =d@  and  cos (dd)  =1. 
Then,  Atf  -fa+u.)de+  dud&-t-dvjj  +ldu-vd9- dvd&Jj  .  (la-5) 

Neglecting  infinitesimals  of  higher  order, 

A,  P,  ~ [(Z+u)d&+dv]T  +[du  +  vdO)  J  (la-6) 

and  \A,?d=L(r*+u)Xd^+zfo+<i)vddZ+CL*'dO l-Zu.vdOi+v\/0J  .  (la-7) 


The  strain,  ,  on  the  mid-thickness  curve  is  defined 
as  the  change  in  length  per  unit  initial  length  of  the  curve 
as  the  tube  wall  deflects.  For  small  angles  the  arc  length 
and  the  cord  length  can  be  taken  as  equal.  Therefore,  the 


strain  is 


ee,/ 


_  M,?f  -  rs  c/<9 
r„d& 


(la-8) 


Substituting  the  value  of  lA,Vl  yields 


_  IC.CHjfdd’ tzCK+U.)vdd^+OL'c/&X-ZU-Vd&i'+  vV^J  -  ^  J-& 
*>'  r9  dO 


(la-9) 


Again  neglecting  infinitesimals  of  higher  order,  the  strain 
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becomes  r<>  [/  +  ~ (u.+  v) ]  e  —  f0 


(la-10) 


O 


(la-11) 


Therefore , 


-  VkL  v 
'&>!  K 


(la-12) 


on  the  mid-thickness  curve. 

The  second  component  of  strain  is  caused  by  the  change 
in  curvature  of  the  tube  wall*  Because  the  tube  is  "thin" 
the  theory  for  bending  of  straight  beams  can  be  used*  The 
strain  due  to  bending  is 

F  &A  ,  (la-13 

where  r  =  distance  from  mid-thickness  and 

=  change  in  curvature  of  tube  wall  mid-thickness 


curve  o 


To  find  the  change  in  curvature  consider  a  point  P 
which  after  a  change  in  temperature  is  at  a  new  location 
P,  as  shown  in  Figure  6* 
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y 


Figure  6.  Displacement  of  a  Point  on  Mid-thickness 
Curve . 

The  point  is  located  by  the  vector  R=  (C>+CL)er+  V&e  (la-14) 

where  &r  and  are  unit  vectors  in  the  radial  and  tan¬ 
gential  directions  respectively.  Taking  derivatives; 

#  «  uer+  ( r0+u)e6+  veff  -  ver 

=  (u-v)er +(rz  +  u+  v)ea  (la-15) 

and  S  -  (a-v)er  +(ul-  V)e&  +  (u+ V)eB  -  fa+u.  4 1 >)  er 

=  (u.-Zv-CL-r0)er+(V+Zu~v)ee  .  (la-16) 
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The  curvature  after  deflection  isfeJ 


i/_  (u-v)(u+2cl-v)-(cI  -Zv-U‘-KXr*  +  u.+  v)  M  1rn 

•  (  71 

Noting  that  re  greatly  exceeds  u,  v  and  their  derivatives, 
we  get  the  approximation 

y  1 

=  r0z-(oL+&)r0  .  (la-18) 

r3 

fo 


The  change  in  curvature  is  the  final  curvature  minus  the 
initial  curvature.  The  tube  is  initially  round  so  the 
initial  curvature  is  l/r0  .  Thus, 


,K=  r/-(u+“)  -  1  =  - 


r* 


6 


(la-19) 


The  strain  due  to  bending  is  then 


(g±u)r 

Cz 


(la-20) 


The  tangential  strain  is  the  sum  of  the  two  components 
developed  above : 


•  -  v  _  (udi  u)  r 

*  C  £* 


(1) 
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Solving  Equation  2a-4  for  a we  have 


a.  -  0^)8o  -  L 1  ( 2a-5) 


and  <Z-  0  +  y)£0  ~~ 

r0 

7. 

£  =  /  . ,  ^  ■*. 
where  /  /<r«  . 

(12) 

Equation  6. 

£<%ydA=0 

(6) 

y  =  (r0+r)cos  9 

and 

dA=tr0+r)  drdO 

(S) 

Substituting  Equations  8  and  2a-l  into  6  and  integrating 
over  the  same  interval  as  before: 

— f(rz+rf  drj[a.cos0+  b(r0+  ^cos^O +VZ.T  [Li+z(n',)fao+  * 
l-u  £  '*-» 

~C)-hi02iSjzo^C^O)cosO]clO  ~  0 ,  (2a-6) 

>1*0  ^ 


r  tt  :  * 

Note  that  J  cosC*e) cos(*9)d0 ^  M  . 


(2a-7) 


Henceforth  in  the  analysis,  terras  which  vanish  in  this  way 
will  not  be  written.  Continuing  the  integration, 

32 

I  [b (r0+rf +y(u,+  V)(z±jj2 - (/+»)&,  (r0+r?’]cir  =  0  (2a-8) 

'  4  * 
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Let  <f-  / 


=  /  + 


tz  c*  > 


a  rs  °-h<*  ~ 

then  b  =  ~A~  [0+v)£,  -  Vjj^+vj  ] 

<*  £  ra 


Equation  15 . 

zo-pv = f  ~f\+el  -&+»)£&)  ^ 

-2-(]+v)8(d)€9  +  zO+v)e(efl  JA 


(13) 


(15) 


The  integral  is  the  sum  of  the  following  terms: 

I  / e*  JA  =  fr.-tr)jrJ[a}+  bVr„+r/^s^j 

^  -to 

x 

=  fWc)  fe71" &+  if  bx(ra+  r)^]d  r 

'i 

=  £~TT GL  lrQ{+  7 TbV»3 ^  f/ * ] 

*  zira^ni  -hir^-tlz  £ 

A  *7T 

sr  J^^r)drfjzJ2.{[i  4~E/»x'!)]  «£  +Znd+'6’>-d]un^  +  *T  t&JcAStfisWd 


IZZT 


'4 


o  °  o 


r% 

= J  (&*)?*  fa  -£)£+  Z  +2"C'+  %«'$** 


-K 


zzi  ]z £  (K+^^fa^+ZnP^ U-H^-i-  tfXVJ 

C  I 
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I 


where  £  =  I+/Ti;z^~^  a>act  /°r,=  ,+  • 

k  ri7r 

€&JA  =  C(rc+>OcIrJ[£i-£)a.<jL0  +  (ra-br)(ul+\/))bcos'Ld}  oS 


/a. 


-JtL  J2Ve. 

•  v  r  J 

=  SJjJT  J (tr+r)l20+  gJaUt  +  CKt+rXUf+vJbJdr 
IZ 

-  ZvillJ  IzclH*  £  +  (u,+  vi)b{$ 

^  Jit 

JV'  f-z{/-hv)£(d)eiclA=  -2(n-y)f  (n,+r)JrJla.£()'h  b(K+r)8lcosz0]d 6 

\ 


=  -2  Oi-y)yr  flz  a.^+  6  ^  r 

4 

=  -ZO^TTrzt  [za.£0  +  b  Toft  £,] 


*/z  ** 

C-zO+tJgfc)  e9  dA = 

/"&  « 

*0  -it/  * 


=  z^jrsi  [zfu.£'  +f(AHibn  <)£„] 

*A  /  ' 


£  >7T 

W  J  Z(f+»j£(a)ZdA  =  2  0+y)f.fc+r)di/y  8^  cos  *6*©)  d-9 

A  ~o 


=  zO-n>)-*Q  t  [2 e*  +Z  <£*  j 
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Summing  all  terms: 


ZO-W- 

7Tg  i£ 


=Zaz+  +  ^,[2  £u?  (^u*  +2  npn  Vn+n'K\ ? 

+Zr[2GUo  f+Cu?  v)bfi]-2(i+»)L&-£a+  b/3r06,] 


-z^[z$Uv£0+^(fiu.»+nV„)Sn  J 


+zO+v)[2£o  +1^£«] . 


(16) 
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APPENDIX  III 


Example 

•  The  following  example  will  illustrate  the  method  of 
solution  and  order  of  magnitude  of  results.  A  302  stain¬ 
less  steel  pipe  10  inches  in  diameter,  with  a  wall  thick¬ 
ness  of  0.5  inches  is  subjected  to  the  temperature 
distribution  shown  in  Figure  3.  Assuming  the  pipe  is  free 
from  restraint,  what  is  the  stress  pattern  and  how  much 
will  the  pipe  bow? 


e,  RADIANS 


Figure  3.  Temperature  of  Pipe. 
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The  first  step  is  to  determine  the  thermal  strain  for 
each  temperature .  The  base  temperature  in  this  case  is 
taken  as  6£°F  because  the  data,  see  Figure  9,  is  given  in 
that  form.  A  harmonic  analysis  is  then  performed  to  deter¬ 
mine  £o >  S,  j  and  S,  .  The  calculations  and  results  are 
given  in  Table  II  and  paragraphs  below. 


TEMPERATURE ,  °F 

Figure  9.  Thermal  Expansion  of  AISI  302  Stainless  Steel 

From  T.  F.  Durham,  R.  M.  McClintock,  and  R.  P.  Reed, 
"Cryogenic  Materials  Data  Handbook",  National 
Bureau  of  Standards,  PB  Report  171$09,  1961. 


TABLE  II 


Harmonic  Analysis 


8 

T(0) 

SCO) 

CO 

0 

0 

§ 

£6^3  sin 

radians 

°F 

X/0* 

x/o ** 

x/or 

0.0 

50 

-20 

-20.0 

0.0 

0.1 

43 

-21 

-20.0 

-6.5 

0.2 

43 

-26 

-21.1 

-15.3 

0.3 

35 

*35 

-20.6 

-23.3 

0.4 

23 

-45 

-14.2 

-43.3 

0.5 

-10 

-62 

0.0 

-62.0 

0.6 

-20 

-30 

24.7 

-76.2 

0.7 

-70 

-125 

73.5 

-101. 4 

0.3 

-250 

-260 

211.0 

*153.0 

0.9 

-295 

-235 

272.0 

-33.2 

1.0 

-300 

-233 

233.0 

0.0 

1.1 

-295 

>235 

272.0 

33.2 

1.2 

-250 

-260 

211.0 

153.0 

1.3 

-115 

-162 

95.3 

131.2 

1.4 

-100 

-150 

46.3 

142.5 

1.5 

-35 

-140 

0.0 

140.0 

1.6 

-75 

-130 

-40.2 

123.3 

1.7 

-50 

-110 

-64.7 

39.2 

1.3 

15 

-50 

-40.5 

29.4 

1.9 

40 

-30 

-23.5 

9.3 

Total 

-2566 

1224.0 

323.0 

!.9V 

E0  =  ^  II  €(e)  =  -  IZ8- 3  Xfo' 

/•d'tf 

E,  «  TslE(e)cose  =  izZA  x  /a 

o 

£,  =  7^  <? (&Js//7  0=?32.3x  /o' 


The  bowing  of  the  pipe  is 


—  Z.^3x/0  /In 


n> 


at 


O IP  —  /a.*-'  S,  -  i&n1  32A  =.  /4.  8  . 


IZ2-4 
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The  tangential  stress  is  given  by 


$e  ^  e  sm  9  ^  > 

therefore,  on  the  mid-thickness  curve,  r=0,  <7g- 0. 

Letting  E=30xl06psi  and  >'=0.296, 

vA.  _  o.z9i(3o*)  _  q  7^6 

/-^  ,-o.o&iz  ~  y'/6  eSj-  • 

At  the  outside  diameter  r=t/2=0.25  inches. 

0 ~Q  -  9.7  cos#  +  32-3  s/v  0  )t 

The  results  are  given  in  Table  III. 

The  axial  stress  is  given  by 

-E{s(&) -]£*+ 0+ +  €,  sw9) ]} 

=  -[-iz&f+fi  +rtfrP))(l2Z'fcos  d  +3Z-ifsme)  ]}  # 

On  the  mid-thickness  curve  r=0  and  [  I  ^  =1» 

At  the  outside  diameter  r=0„25  and  [/+  3  -1-055- 

Axial  stress  was  calculated  for  each  position  and  the 
results  are  given  in  Table  III  and  also  shown  in  Figure  10. 
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TABLE  III 


■ 


Tangential  and  Axial  Stress 


6 

,psi 

, psi 

0i  , psi 

radians 

at  r=0.25 

at  r=0 

at  r=0.25 

0.0 

596 

4,200 

5,700 

0.2 

575 

4,800 

8,700 

0.4 

333 

“4,200 

-3,000 

0.6 

-23 

”16,500 

-16,800 

0.8 

“390 

14,400 

13,500 

1.0 

-596 

11,400 

9,300 

1.2 

”575 

4,200 

2,400 

1.4 

”333 

-13,800 

-15,000 

1.6 

23 

2,700 

3,000 

1.8 

_ m _ 

600 

_ 2,700 _ 

Figure  10.  Stress  Pattern  of  Pipe. 
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The  result  of  the  first  assumption,  that  <Tr  and  Tr$ 
are  negligible,  is  that  equilibrium  cannot  be  satisfied. 
The  values  of  these  stresses  are  now  estimated  using  the 
conditions  of  the  example.  A  section  is  cut  out  as  shown 
in  the  free  body  diagram,  Figure  11.  These  are  the  posi¬ 
tions  of  maximum  and  zero  tangential  stress. 


y 


Figure  11.  Shear  Stress 

The  stress  distribution,  ,  causes  a  net  counter-clock¬ 

wise  moment  of  25.7  in-lbs.  To  restore  equilibrium, 
stresses  0^  and  7?re  are  required.  Summing  moments  about 
the  point  P, 

7ro*Mih  ■  /0-2  Ps‘  ■ 

The  tangential  stress,  <rQ'  ,  is  also  equal  to  10.3  psi. 

An  indication  of  the  magnitude  of  the  radial  stress 
required  to  maintain  equilibrium  can  be  obtained  by  inves¬ 
tigating  the  condition  at  0=195°,  the  cold  side  of  the  pipe. 
There  the  axial  stress  is  about  10,000  psi  and  will  be 
taken  as  constant.  Figure  12  shows  the  action  of  axial  and 
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radial  stress. 


<%=  10,000  pS( 


Figure  12.  Radial  Stress  Component  Caused  by 
Axial  Stress. 


The  radius  of  curvature,  R,  was  obtained  from  the 
bowing.  Summing  forces  in  the  radial  direction, 

‘R(Trd6- {(T^dd  ~0  . 

Therefore,  <V '  ^§'^0°°°  ^  =  1.27  psi’  the  radial 

stress  caused  by  axial  stress. 

The  maximum  radial  stress  component  resulting  from 
tangential  stress  can  be  estimated  by  considering  the  stress 
condition  at  S- dp  as  shown  in  Figure  13. 
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Figure  13 .  Radial  Stress  Caused  by  Tangential 
Stress. 

Taking  the  sum  of  the  forces  in  the  radial  direction 
we  get, 

z<rr  r0J&- a  ^C\<rd  +  %' )6d 
^ r  ~  °i/§j  +  /0)  ~  776  Psc  • 

The  total  radial  stress  is  the  sum  of  the  component 
caused  by  axial  stress  and  the  component  caused  by  tangen¬ 
tial  stress  (  G~r  ~  •=  9  ps  <  ) . 

It  is  seen  from  the  above  calculations  that  op.  and 

are  small  when  compared  with  other  stresses,  and  there¬ 
fore,  neglecting  them  is  realistic. 
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